Noncoplanar magnetic ordering driven by itinerant electrons on the pyrochlore lattice 
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Exchange interaction tends to favor collinear or coplanar magnetic orders in rotationally invariant 
spin systems. Indeed, such magnetic structures are usually selected by thermal or quantum fluc- 
tuations in highly frustrated magnets. Here we show that a complex noncoplanar magnetic order 
with a quadrupled unit cell is stabilized by itinerant electrons on the pyrochlore lattice. Specifically 
we consider the Kondo-lattice and Hubbard models at quarter filling. The electron Fermi 'surface' 
at this filling factor is topologically equivalent to three intersecting Fermi circles. Perfect nesting 
of the Fermi lines leads to magnetic ordering with multiple wavevectors and a definite handedness. 
The chiral order might persist without magnetic order in a chiral spin liquid at finite temperatures. 



Magnets with geometrical frustration have fascinated 
physicists for more than a decade as models of strongly 
interacting systems. A prominent example of magnetic 
frustration in high dimensions is the nearest-neighbor an- 
tiferromagnet on the pyrochlore lattice shown in Fig.[lja). 
For classical Heisenberg spins, strong geometrical frustra- 
tion prevents the magnet from settling in a long-range 
magnetic order even at zero temperature [l|. Instead, 
the ground state retains a finite zero-point entropy and 
is susceptible to small perturbations such as anisotropics, 
further-neighbor couplings, and dipolar interactions. In 
real compounds, magnetic frustration is often relieved 
when spins couple to other degrees of freedom. An inten- 
sively studied case is spin-lattice coupling in chromium 
spinels where the magnetic transition is accompanied by 
a lattice distortion [2-|j|. Another well known example 
is magnetic phase transition faciliated by orbital order in 
frustrated spin-orbital systems ^5-8]. 

Metallic pyrochlore magnets such as Kondo-lattice or 
double-exchange models pose a different challenge for 
theorists due to the nonlocal nature of the electron- 
mediated interactions. Indeed, despite considerable effort 
[9l4l2|. a complete picture of the phase diagram is still 
lacking. An early study of the double-exchange model 
on pyrochlore lattice revealed a rich mean-field phase di- 
agram [9| . However, the calculation only considered mag- 
netic orders which preserve the lattice translational sym- 
metry, hence precluding magnetic structures with mul- 
tiple wavevectors. A recent Monte Carlo study of the 
same model provided an unbiased phase diagram at large 
Hund's coupling [llj. In particular, a peculiar state with 
electronic phase separation was observed. Magnetic or- 
ders in the weak-coupling regime remain unclear. 

Magnetic ordering in frustrated metallic systems de- 
pends in an intricate way on the underlying electron 
Fermi surface. At small filling factors, localized spins 
interact with each other through a long-range oscilla- 
tory Ruderman-Kittel-Kasuya-Yosida (RKKY) interac- 
tion mediated by the electrons. A Monte Carlo sim- 
ulation of pyrochlore spin-ice with RKKY interaction 
showed that the sign of the effective Curie- Weiss con- 
stant as well as that of nearest-neighbor coupling vary 



with the electron Fermi wavevector 10]. This in turn 
determines magnetic ordering at low temperatures: the 
magnetic state of individual tetrahedra evolves from the 
all-in all-out structure to the 2-in 2-out ferromagnetic 
state with increasing electron density. 

The geometry of the electron Fermi surface also plays 
an important role in determining the magnetic instability 
of itinerant systems, particularly at commensurate filling 
factors. A canonical example is Neel ordering caused by 
perfect Fermi surface nesting that occurs at a half-filled 
square-lattice Hubbard model. While the resulting spin 
structure is collinear in bipartite square lattice, it was 
shown that nesting effect on triangular lattice leads to 
a rare noncoplanar magnetic order at filling factors 3/4 
and 1/4 [13|-[l5|. This chiral magnetic structure induces 
a spontaneous quantum Hall effect in the absence of ex- 
ternal magnetic field. Although a long-range magnetic 
order cannot survive thermal fluctuations in two spatial 
dimensions, the discrete chirality order persists up to a 
finite temperature [13|. 

In this paper we demonstrate that the ground state 
of an isotropic Kondo-lattice model on the pyrochlore 
lattice is magnetically ordered with a noncoplanar spin 
structure and a definite chirality. The conclusion also 
applies to Hubbard model on the pyrochlore lattice at 
the mean-field level. The noncoplanar magnetic order 
stems from a weak-coupling instability caused by perfect 
nesting of Fermi 'circles' at quarter filling; the quadru- 
pled magnetic unit cell contains 16 spins. In contrast to 
the triangular-lattice model |13j . this noncoplanar spin 
order does not support a spontaneous Hall insulator be- 
cause of a trivial Berry phase acquired by the electrons 
when traversing a tetrahedron. The magnetic structure 
itself, on the other and, has a definite handedness and is 
characterized by a nonzero chiral order parameter. 

It should be pointed out that although chiral magnetic 
orders have been reported in metallic pyrochlore oxides 
such as Pr2lr207 [l6|, the non-coplanarity of spins in 
these so called spin-ice compounds are mainly created by 
a strong uniaxial anisotropy. The electron-driven mag- 
netic instability in pyrochlore spin-ice at quarter filling 
will be briefly discussed at the end of the paper. 



We begin with the isotropic Kondo-lattice Hamiltonian 
on the pyrochlore lattice: 
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The first term describes electron hopping between 
nearest-neighbor sites, t is the hopping integral, and c|q, 
creates an electron with spin a on site i. The second 
term represents the superexchange interaction between 
neighboring localized spins S^. The itinerant electrons 
interact with the localized spins through an on-site ex- 
change coupling K\ the vector (Tai3 denotes the three 
Pauli matrices. We consider the classical limit of local- 
ized spins with magnitude |Si| = 5* 3> 1. In this classical 
limit, the electron spectrum is independent of the sign of 
K and the eigenstates of opposite signs are connected by 
a global gauge transformation [13]. 

We first consider the tight-binding spectrum of elec- 
trons on the pyrochlore lattice. In momentum space, the 
hopping Hamiltonian [first term in Eq. ([1])] is expressed 



Ht = Em.nEk,a^m"(k)cJ„a(k)c„„(k), whcre m, n 
sublattice indices and the hopping matrix is 


T(k) =-2t 


1 Cy^ Czx Cxy \ 

Cyz U Cxy Czx 

Czx Cxy U Cyz 
\ Cxy Czx Cyz U / 


(2) 



Here we have introduced Cab — cos [{ka + fcf,)/4] and 
Cab = cos [{ka — kb) / 4] and set the length of a cubic unit 
cell a = 1. This hopping matrix can be exactly diago- 
nalized, giving rise to a spectrum shown in Fig. [Ijd). It 
consists of two degenerate flat bands e^*** — 2t and two 
dispersive bands: 



2t ITVI + Qk 



(3) 



where Qk = cos ^ cos ^ -I- cos ^ cos ^ -I- cos ^ cos ^. 
These two branches touch each other along the diagonal 
lines of the square faces of the Brillouin zone [Fig. IHb)]. 

At quarter filling the lowest band e^ is completely 
filled. The corresponding Fermi 'surface' is determined 
by e^ = e^ and consists of diagonal lines 27r(a::, 0, 1), and 
so on, on the square faces of the zone boundary. The 
Fermi energy lies exactly at Ep — —2t. Noting that 
27r(a;, 0,1) = 27r(l — 2^,1,0) module a reciprocal lattice 
vector, these diagonal lines are topologically equivalent 
to three closed loops, or circles, intersecting with each 
other at the X points of the Brillouin zone [Fig. [ijb-c)]. 
In the vicinity of these Fermi lines, the spectrum has the 
form 6]^ « Ep:kCip x |k_L|, where kj^ is the perpendicular 
component of electron momentum. 

The reduction of Fermi surface to Fermi lines at quar- 
ter filling also makes it possible to realize perfect nesting 
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FIG. 1: (a) A conventional cubic unit cell of the pyrochlore 
lattice. I and II denote the two types of tetrahedra with 
different orientations, (b) First Brillouin zone of the fee lat- 
tice. At quarter filling, the electron Fermi 'surface' consists 
of diagonal lines on the square face of the Brillouin zone, (c) 
Topologically, these Fermi lines are equivalent to three closed 
loops interceting at the X points of the Brillouin zone, (d) 
Band structure of the tight-binding model on pyrochlore lat- 
tice. The segment XW corresponds to one-quarter of a Fermi 
line at quarter filling. 



with a few wavevectors. Indeed, the three Fermi circles 
can be perfectly nested by wavevectors Qi — 27r(l,0,0), 
Q2 = 27r(0,l,0), and Q3 = 27r(0,0, 1), respectively. For 
example, states on the two Fermi lines 27r(l,0, z) and 
27r(0, l,z), both belonging to the blue Fermi circle in 
Fig-[Uc), differ by Q3 up to a reciprocal lattice vector. 
The Van Hove singularity of the density of states at the 
Fermi level gives rise to a logarithmically divergent Lind- 
hard susceptibility: x(q) ~ — log|q— Q^|, indicating a 
weak-coupling instability caused by the nesting effect. 

To investigate possible magnetic ordering caused by 
perfect nesting of the Fermi lines, we consider spin orders 
which can be expanded as 



S, = S™(r) =^M„,^, 

r;=0 



a'i-Q.ri- 
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where m is the sublattice index and r denotes the posi- 
tion of the Bravais (fee) lattice point. In order to accom- 
modate, e.g. ferromagnetic ordering, we have included a 
Qo = component in the above expansion. The unit cell 
of this multiple- Q magnetic order is extended to the con- 
ventional cubic unit cell which contains 16 inequivalent 
spins. 

Next we include the on-site Hund's exchange which 
couples electron with momentum k to that with k-|- Q^. 




FIG. 2: (a) A unit cell of the noncoplanar magnetic order 
viewed from the c axis. Here we choose fii = a, n2 = b, and 
fis = c in Eq. (0 . (b) Chiral structure of the same magnetic 
order projected onto a (111) kagome plane. The and ® 
symbols denote spin component coming out of and into the 
plane, respectively. Same handedness (right-handed in this 
case) is observed in other kagome planes. 



After Fourier transformation, the on-site coupling term 
has the form 



^K-EEE'cL(k+Q5)x 



(m) 



(k+Q^). (5) 
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Here the prime indicates that the momentum summa- 
tion is restricted to the reduced Brillouin zone, and the 
coupling coefficients are given by 



M 



(m) 



= K A^^, (; [MjrtX ■ O-afl) 6 



^(Qe-Q,)-d„ 



(6) 



where d„i denotes basis vectors of the pyrochlore lattice. 
Momentum conservation of the coupling is encoded in 
the symbol Aj,,_^ which is symmetric with respect to 
the first two indices. The nonzero components include 

^01,1 — Ao2,2 — Ao3,3 = A234 ^ ^31,2 = Ai2^3 = 1. 

For a given set of Fourier components M™ ,,, diago- 
nalization of the electron Hamiltonian Ht + Hk gives a 
total of 32 energy bands. Due to perfect nesting of Fermi 
lines, the first 8 bands are expected to be separated from 
other high energy states by an energy gap. At quarter 
filling, the ground-state energy of electrons E^ is obtained 
by filling the lowest eight bands. However, the determi- 
nation of the minimum- energy state is not straightfor- 
ward because of the large number of variables required 
to describe the spin structure. Indeed, even restricting 
to magnetic orders which can be expressed in the form 
of Eq. (JH), one still needs 32 variables to specify the 16 
classical spins in a quadrupled unit cell. Here we employ 
the simulated annealing algorithm to minimize the en- 
ergy functional Ee (Mm,,,). For Hund's coupling below 
a critical value Kc ~ 4.32 1 and J = 0, we obtained the 
same noncoplanar magnetic order shown in Fig. [2] using 
different initial configurations and rates of decreasing the 
effective temperature. For K > Kc ferromagnetic order 
with Si = S h takes over and becomes the ground state. 
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FIG. 3: (a) Electron band structure corresponding to the non- 
coplanar magnetic order shown in Fig. (2] The Hund's cou- 
pling is i^ = f and superexchange constant J — 0. (b) The 
reduced Brillouin zone with high-symmetry points and lines. 



To characterize the noncoplanar spin structure, we in- 
troduce three order parameters Li = (So + Si — S2 — 
83)745', and so on [4|, which measure the staggered mag- 
netizations of a tetrahedron; the spin subscript indicates 
the corresponding sublattice. A unit cell of the magnetic 
order is shown in Fig. [2ja), where the magnetic states of 
type-I tetrahedra are described by order parameters 



Li 



V5 ' 









(7) 



Here n^'s are three mutually orthogonal unit vectors. 
Remarkably, the noncoplanar magnetic order is a simul- 
taneous ground state of the inter-site exchange interac- 
tion. To see this, we note that the superexchange term in 
Eq. dl]) can be recast into Hj = (J/2) J^s |M^| , where 
M^ = J2m=o ^™ ^^ *^® total magnetization of a tetrahe- 
dron. For the magnetic order shown in Fig. [5J the total 
spin of every tetrahedron, including both types, vanishes 
identically, hence minimizing the exchange energy. This 
result also indicates that the ferromagnetic order favored 
by a large Hund's coupling K is suppressed by the anti- 
ferromagnetic inter-site exchange. 

Fig. EJa) shows the electron band structure corre- 
sponding to the noncoplanar magnetic order. The spec- 
trum is at least doubly degenerate thanks to a combined 
lattice translation and spin rotation symmetry. Due to 
the nesting effect, a charge gap opens between the lower 
four pairs of bands and other high energy branches. It is 
known that the electron experiences a fictitious magnetic 
field proportional to the spin chirality Xijk = Si • Sj x Sfc 
[l6| when hopping around a triangular loop. Although 
Xijk is nonzero in the noncoplanar structure, the contri- 
butions from the four triangular faces of a tetrahedron 
cancel each other due to a monopole-like spin configu- 
ration. The noncoplanar magnetic order thus does not 
exhibit a spontaneous quantum Hall effect. 

The magnetic structure itself, nonetheless, breaks the 
chiral symmetry. This is best illustrated by projecting 
the triple- Q magnetic order onto a (111) kagome plane 
which shows hexagons and triangles with a definite hand- 



edness [Fig.[2jb)]. More specifically, we introduce a scalar 
chirality for a tetrahedron in the antiferromagnetic state: 



XKI 



:= 3\/3 Li • La X L3 



(8) 



The noncoplanar magnetic order shown in Fig. [5] has 
the maximum chirality (xk) = +1, where the average 
is over both types of tetrahedra. This is to be contrasted 
with the so called all-in all-out structure [Fig. SJ^a)] in 
which the two types of tetrahedra have opposite chirali- 
ties (xi) — — (xii); the average chirality is thus zero. 

The model described by Eq. (U]) with a negative K was 
recently investigated in Ref. 11 using Monte Carlo sim- 
ulations. The authors also considered the special case 
of quarter filling. However, since the simulations were 
performed assuming a large Hund's coupling, the non- 
coplanar magnetic order shown in Fig. [5] was not ob- 
served in their simulations. Although the fine structure 
of the electron Fermi surface is absent in small finite sys- 
tems, the noncoplanar magnetic order could be found in 
Monte Carlo simulations at small Hund's coupling using 
techniques such as twisted boundary conditions [11,]. 

We now turn to the spin-ice compounds where the ro- 
tational symmetry is explicitly broken by a strong easy- 
axis anisotropy Ho — —D J2ii^i'^m)'^- The effective de- 
grees of freedom become discrete Ising variables as spins 
are forced to point along the local (111) directions. We 
start with the J = limit. At filling factor ly = 1/4, 
simulated annealing yields a ground state with the all-in 
all-out magnetic order shown in Fig. UJ^a) . The fact that 
every tetrahedron has total spin Mg| — makes the all- 
in all-out structure a simultaneous ground state of the 
antiferromagnetic inter-site couplings. 

In spin-ice compounds, however, the nearest-neighbor 
exchange has a ferromagnetic sign, i.e. J = —\J\- As- 
suming D ^ \J\ ^ t, K, numerical minimization gives 
a single-Q 2-in 2-out structure [Fig. |4ljb)] . The magne- 
tizations of type-I and type-II tetrahedra point along a 
and b axes, respectively, with their sign alternating be- 
tween adjacent ab planes. There is a total of six degen- 
erate single-Q states related by time reversal and lattice 
translation symmetries. 

Both magnetic orders shown in Fig. 2] were obtained in 
a recent Monte Carlo study on a metallic pyrochlore spin 
ice [1^ . Instead of tackling the Kondo-lattice model di- 
rectly, the simulations were performed assuming a long- 
range RKKY interaction between localized spins. Ne- 
glecting the superexchange J, their simulations yielded 
a 2-in 2-out structure at quarter filling [lO|, which is in 
disagreement with our result (the all-in all-out state). 
This discrepancy could be attributed to the fact that the 
RKKY interaction is derived by integrating out the elec- 
trons with a spherical Fermi surface. The approach thus 
neglects the subtle effect of Fermi surface nesting. 

In summary, we have studied a model of itinerant elec- 
trons interacting with localized classical spins on the py- 
rochlore lattice. At quarter filling, we showed that a 




FIG. 4: (a) All-in all-out structure, (b) 2-in 2-out magnetic 
order with wavevecotr Q3 = 27r(0, 0, 1). The and ® sym- 
bols again denote spin component coming out of and into the 
plane, respectively. 



noncoplanar magnetic order emerges as a ground state 
of the rotationally invariant Hamiltonian. The magnetic 
structure characerized by a nozero order parameter XB 
also breaks the chiral symmetry. Since the chiral transi- 
tion needs not coincide with the magnetic ordering from 
the symmetry viewpoint, a chiral phase with (xki) 7^ 
but no spin order (L.^) = could occur at finite tem- 
peratures. A similar case was recently observed in spin 
ice compound PralraOr |17|. Future studies, especially 
large-scale Monte Carlo simulations, could demonstrate 
the noncoplanar magnetic structure numerically and ex- 
plore the possibility of the partially ordered chiral phase. 
Finally, we showed that when the localized spins are sub- 
ject to strong uniaxial anisotropy, the on-site Hund's 
coupling favors a uniform all-in all-out spin structure, 
whereas the single-Q 2-in 2-out magnetic order is stabi- 
lized by the inter-site ferromagnetic exchange interaction. 
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